A droplet near the critical point: the divergence of Tolman's length 
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Abstract 

Application of "complete scaling" [Kim et al, Phys. Rev. E 67, 061506 
(2003); Anisimov and Wang, Phys. Rev. Lett. 97, 25703 (2006)] to the 
interfacial behavior of fluids shows that Tolman's length, a curvature cor- 
rection to the surface tension, diverges at the critical point of fluids much 
more strongly than is commonly believed. The amplitude of the diver- 
gence depends on the degree of asymmetry in fluid phase coexistence. A 
new universal amplitude ratio, which involves this asymmetry, is intro- 
duced. In highly asymmetric fluids and fluid mixtures the Tolman length 
may become large enough near criticality to be detected in precise exper- 
iments with microcapillaries and in simulations. 
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For decades the behavior of Tolman's length has remained one of the most 
controversial issues in mesoscopic thermodynamics. The Tolman length (5) is 
defined as a curvature-correction coefficient in the surface tension (cr) of a liquid 
or vapor droplet pQ: 

a(R)=a 00 (l-^ + ..)j, (1) 

where R is the droplet radius, taken equal to the radius of the surface of tension 
E], and Coo is the surface tension for the planar interface. The physical 
origin of the curvature correction is a difference between the equimolar surface 
and the surface of tension; the difference is phenomenologically associated with 
asymmetry in fluid phase coexistence. In symmetric systems, such as the lattice- 
gas model and the regular-solution model, the difference between the equimolar 
surface and the surface of tension vanishes and the Tolman length does not 
exist [310]. The curvature dependence of the surface tension is important for the 
description of nucleation phenomena and fluids in microcapillaries or nanopores. 
However, while square-gradient theories gave consistent results in mean- field 
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approximation 01 El El j the actual critical behavior of Tolman's length is not 
certain, either in sign or in behavior [7| E] . Fisher and Wortis [2| and Rowlinson 
jS] have predicted a very weak algebraic divergence at the critical point with 
an exponent —0.065. Since this prediction is supported by an exact expression 
|H] obtained for the Widom-Rowlinson "penetrable-sphere model" UJ, it has 
become commonly accepted. 

In this Letter, I show that the application of "complete scaling" |10l 
which properly describes vapor-liquid asymmetry, to the problem of Tolman's 
length yields a much stronger algebraic divergence at the critical point, with 
an exponent —0.304. The amplitude of this divergence depends on the degree 
of asymmetry in fluid phase behavior. In highly asymmetric fluids and fluid 
mixtures this divergence may be strong enough to be detected in accurate ex- 
periments and simulations. 

The generalized Laplace- Tolman equation reads E| • 

2a 00 ( 25 \ 

Pa -P,=-^(i-- + ...y (2 ) 

where P a — Pp is the difference between the pressure inside the droplet (phase 
a) and outside (phase (5). The difference P a — Pp in Eq. (2) can be separated 
into two parts, namely, a symmetric part, which does not depend on which fluid 
phase is inside or outside the droplet, and an asymmetric part, i.e., P a — Pp = 
(P a — Pp) sym + (P a — -P/3) asym - Since Tolman's length is associated with the 
asymmetric part only, one can conclude that 2a(oo)/R = (P a — Pp) s m , and 
thus in first approximation obtain the expression 

In scaling theory [2] the surface tension in a three-dimensional system can 
be estimated as a product of the " critical part" AP cr of the grand potential per 
unit volume (fl/V = —P), which scales as AP cr ~ £~ 3 , and the characteristic 
"thickness" of the interface, the correlation length of density fluctuations £. The 
correlation length diverges asymptotically as £ w £J |AT| " with the universal 
critical exponent v ~ 0.630 and the amplitude (below the critical point) 

ocT 2 



AT 



while the surface tension vanishes at the critical point as Coo w ctq 

|12| . Here and below, the circumflex indicates dimensionless variables, such 
as AT = (T-T c )/T c , AP = (P - P c ) / Pc k B T c , A/t = (n- Pc )/k B T c , and 
A/? = (p — p c ) I p c , with /cb being Boltzmann's constant, while the subscript 
c indicates the critical value of temperature, pressure, chemical potential, and 
molecular density, correspondingly. 

The function AP cr can also be separated into symmetric and asymmet- 
ric parts. For a planar interface, Coo/ (AP cr ) sym £ w c CT , where c a ~ 3.79 
is a universal constant. The mean-field value of c a is 16/3 13 . Therefore, 
(Pa - P ) sym R scales as (AP cr ) sym ( and (P a - Pp) asym R 2 scales as (AP cr ) asym f 
Hence, it follows from Eq. (4) that 
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26 _ (APcr) asym 

€ ~ * (AP cr ) sym ' W 

where c$ is another universal constant. As shown below, the mean-field value 
of c«5 is 5/6, while the scaling-theory value is estimated as ~ 0.6 — 0.7. 

Close to the critical point, AP cr can be approximated as AP cr w x(A/i) 2 , 
where the derivative x — {dp/dfi) T = (d 2 P/dp 2 ) _, the susceptibility with re- 
spect to density fluctuations, is taken for liquid or vapor phase (depending on 
which would be considered as a-phase) at the corresponding phase boundary. 
The susceptibility x becomes asymptotically symmetric close to the critical point 
but exhibits growing asymmetry upon departure from the critical point, contain- 
ing two terms, symmetric Xsym (equal for both phases) and asymmetric Xasym 
(having opposite signs for the different phases). Since (AP cr ) m rs x sy m (A/x) 2 
and (AP cr ) asym w Xasym ( A pf E q- (4) becomes 

S Asym 

The vapor-liquid asymmetry in fluid criticality is most appropriately treated 
by so-called "complete scaling", originally introduced by Fisher and Orkoulas 
[Tl] and further elaborated by Kim et al. ^] and by Anisimov and Wang . 
Specifically, it was shown 11 that with an appropriate choice of the critical value 
of the entropy, p c S c — {dP/dT) cxc , where the subscript cxc means vapor-liquid 
coexistence, the asymmetry in fluid criticality is governed in first approximation 
by only two coefficients, a 3 and 62, in the linear mixing of three physical fields 
(pressure, temperature, and chemical potential) into two independent scaling 
fields, the "ordering" field hi and the "thermal" field /12 according to 



hi = A/t + a 3 



(6) 



h 2 = Af + b 2 Afi, (7) 
where the derivative IdP/dTj is taken at the critical point. The corre- 

V / cxc 

sponding field-dependent thermodynamic potential /i 3 (hi, h 2 ) is also a linear 
combination of the physical fields. In traditional ("incomplete") scaling for flu- 
ids, the nontrivial coefficient a 3 , responsible for the pressure mixing into the 
ordering field, is assumed to be zero EH • 

The total susceptibility x = Xsym + Xasym i s ) m general, a combination of 
three symmetric scaling susceptibilities: strong Xi — (d 2 h 3 /dhf), weak X2 — 
(d 2 h^,/dh^, and cross X12 

= (d 2 h 3 /dhidh 2 ) EDUS]. A detailed derivation of 
the susceptibility, based on " complete-scaling" , has been made by Kim et al. 
jl()| . As a lower-order approximation of the result given in ref. jl()| . one can 
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obtain 
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Ap J Xi + + 2(1 + 03)62X1 



(8) 



When a 3 = (traditional "incomplete" scaling), Eq. (8) becomes identical to 
the result obtained in ref. |16| . 

By using expressions for the scaling susceptibilities along the vapor-liquid 
coexistence (hi — 0) obtained in refs. |1UI I16| . one can derive from Eq. (8) 



X 



sym 



AT 



and 







-7 


3a 3 
1 + a 3 




Xasym. ~ "^^O 


AT 




AT 



i + rv 



2b 



AT 



I3B 



AT 



l-a-P 



(9) 
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where T is the asymptotic critical amplitude of the susceptibility, Bq is the 

asymptotic amplitude in the expression Ap sa ±Bq AT , 7 ~ 1.239, (3 ~ 0.326, 

and a ~ 0.109 are universal "Ising" critical exponents; Tj" is the first symmet- 
ric ("Wegner") correction amplitude, and 6 ~ 0.5 is the "Wegner" correction 
exponent; ± corresponds to liquid and vapor densities, respectively. 
By expanding the ratio 25/£ w — c$ \XasYrnlx s one obtains 
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Hence, Tolman's length contains two diverging terms, namely, 
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where =F corresponds to a liquid droplet or a bubble of vapor, respectively. 

The behavior of Tolman's length, given by Eq. (12), differs in an essential 
way from the results obtained by previous investigators E] E| : a new term, 

, emerges from the complete-scaling analysis. The first term in Eq. 
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(14) diverges more strongly, since /3 



-0.304, whereas 1 — a ~ (3 



-0.065. In "incomplete scaling", when a 3 = 0, Tolman's length diverges very 



weakly, namely as 
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0JIHI. However, the "incomplete 



scaling" result may be valid for some specific systems, such as the "penetrable 
sphere model" jSJ |S] , since that model has an exact symmetry axis on which the 
chemical potential is an analytic function of temperature |17j . 
Equation (12) can be further re-expressed as 
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(13) 



with use of the universal ratio c = (3B^/T Q A Q : 
A (pS) = f (CV) cr dT/fc B T a - ^ / (1 - a) 



1.58 ^Hl an d by introducing 

. l-a 
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mensional entropy (per unit volume) deviation from its critical value; Aq is the 
amplitude of the critical (fluctuation-induced) part of the isochoric heat capacity 
(CV) cr in the two-phase region The higher-order term B cr \AT\ , known as 
the "fluctuation- induced critical background", does not contribute to the diver- 
gence of Tolman's length, since 1 — f3 — v ~ 0.044. However, this term is needed 
to satisfy crossover to the mean- field regime ^T]; moreover its contribution is 
of lower order than that from the " Wegner" correction [j3 + 9 — v ~ 0.2). 

In mean-field approximation, where c = 2, f3 — 1/2, a = 0, and v = 1/2, 
Tolman's length remains finite: 
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(14) 



where £ is the mean-field amplitude of the correlation length and c$ is the 
mean- field value of the constant cs; p,^ = i^d l+1 pj 'dp^dT^j are the dcriva- 

— a 1/2 

tives of the dimensionless chemical potential, Bq = (A30/6A11) is the mean- 
field amplitude of the density coexistence curve, and ACV/^b = 3p\i/p3o is 
the mean- field heat-capacity discontinuity {Fij\. The leading asymmetry coeffi- 
cients, which are assumed to be unaffected by fluctuations, can also be expressed 
through the derivatives p,^ 1 1 1 j as 



«3 



l + a 3 



2 A21 

3 An 



1 A40 
5 A30 



&2 = 



1 ( M21 1 A40 

5 A30 



M11 VM11 



Substituting Eqs. (15) into Eq. (14), one obtains the compact result 
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By comparing this with the mean-field result of Fisher and Wortis 0] and 

Giessen et al. 0, namely, <5/£ = ±TjoA 4 o/ 12 A30 = ± ( 6 AnA3o) A4o/ 12 A30> 
one obtains the mean- field value C5 = 5/6. For the (mean- field) van der Waals 
fluid, where 03 = 1/4, 6 2 = 4/45, p 2 i = 0' An = 9/4, A40 = — A30 = — 27/8, and 
Bq = 2 |TI], one obtains (5 = =f£ /6. If one adopts ca = 1, <5 = =f£ /5. Hence, 
Tolman's length in the van der Waals fluid is negative for a drop of liquid and 
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positive for a bubble of vapor, in agreement with the conclusions of Fisher and 
Wortis |U and Blokhuis and Kuipers [B|. The same sign is anticipated for real 
asymmetric fluids when a 3 is large and positive 

It is also interesting that p 21 does not enter in the mean-field result for 
Tolman's length, which depends only on one asymmetry amplitude /i 40 , the 
fifth-order coefficient in the expansion of the Helmholtz energy per unit vol- 
ume in powers of density. When 0,3 = 0, Eq. (15) gives 62 = p. 4o /10fi 30 p, 11 
and the mean-field result for Tolman's length given by Eq. (16) remains un- 
changed. However, in the scaling regime this asymmetry term, renormalizcd 
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and 



by fluctuations, produces two different singularities, namely, oc 

. A renormalization-group treatment of the fifth-order term 
in the initial Hamiltonian, with this fifth-order term regarded as a "non-Ising 
asymmetry", produces a new critical exponent 65 = 1.3 |20| and results in a 
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nondivergent contribution to Tolman's length, oc 
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H. Ac- 



tually, the fifth-order term in the initial Hamiltonian also generates a singularity, 
, making a special treatment of "non-Ising asymmetry" practically 
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irrelevant. 

In order to obtain an estimate for the universal constant cs in the scaling 
regime, one can consider an approximate square-gradient-theory expression for 
Tolman's length E], namely, 



(Xq ~ Xfi) 
1 a - Pfi) 



(17) 



By calculating the difference of the susceptibilities via Eqs. (9) and (10), and 



by using (p a - p g ) « 4B 
a /k B T cPc ^A^ ~ 2.25 [H] 
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and the universal ratios 
0.47. 



Coo ~ Co 

, TqAq/Bq ~ 0.21 18 1, one obtains c s 
Since in mean-field approximation Eq. (17) gives cs = 2/3 instead of the "exact" 
value cs = 5/6 0] E] or instead of (also plausible) cs = 1, the cs value must 
be corrected at least by a factor of 5/4, yielding cs — 0.6 or by 3/2, yielding 
c 5 ~ 0.7. 

One can also notice that Tolman's length can be thermodynamically ex- 
pressed through the amount of excess adsorption, r s , at the surface of ten- 
sion as S = r s / (p a — pp) P E|- Assuming that the excess (asymmetric 
part) density at the surface of tension is associated with the "singular diam- 
eter" (the deviation of the mean density from the critical isochore), p d — 1 = 

{Pa + Pp) / 2 Pc ~ 1 ~ a 3 /(l + as) (Ap) 2 + 6 2 A (pS) [IH, in such a way that 
r s w — (p d — 1), one obtains 5 w — £ (p d — 1) / (p a — pg). If one adopts 
cs = 2/3, this estimate asymptotically agrees with Eq. (13). For the van 
der Waals fluid such an assumption yields <5/£ = ^fD/B Q = =Fl/5j where 



D = (p d -l)/ 



AT 



= — (3/5) PuPm/ = 2/5 is the mean-field rectilinear 
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diameter of the coexistence densities [F3\, in agreement with Eq. (16) if eg = 1. 

The exponent combination (3 — v, which is responsible for the leading di- 
vergence of the Tolman's length, arises in several contexts in surface thermody- 
namics [21] ■ In particular, the total adsorption r ~ £ (p a — pg\ diverges at the 

critical point as AT . Interestingly, in mean-field approximation, where, 
formally, f3 — v = 0, the adsorption, instead of being finite, becomes logarith- 
mically divergent [21 I21| . Whether a similar logarithmic divergence may also 
persist in the mean-field approximation for Tolman's length in the critical region 
remains unclear. 

A straightforward application of the present results to the problem of Tol- 
man's length in asymmetric binary fluids can be made by a generalization of 
"complete scaling" [221 with use of the "isomorphism principle" JHj. In par- 
ticular, one can use the expression for Tolman's length given by Eq. (13) for 
"incompressible" liquid mixtures replacing Ap by Ax = (x — x c ) /x c , the devi- 
ation of the molar fraction x from its critical value x c . 

The divergence of Tolman's length may be large enough to be detected in well 
designed accurate experiments with microcapillaries and in simulations. Predic- 
tions for Tolman's length in several fluids exhibiting strong vapor-liquid asym- 
metry are shown in Fig. 1. In less asymmetric fluids with small-size molecules, 
such as nitrogen or methane, the coefficient 03 is too small |11] to make Tolman's 
length growing to a nanoscale in any reasonable proximity to the critical point. 
However, in highly asymmetric fluids and fluid mixtures, such as high-molecular- 
weight hydrocarbons and their mixtures with small-molccular-volume species, 
or in ionic fluids, the Tolman's length can easily reach a nanoscale. From Eq. 
(13) with eg = 2/3, for the liquid mixture n-hexadecane+nitrobenzene, one ob- 



tains —26 ~ 5 nm at 



AT 



= 1CT 4 , while £ « £0 |ATp 063 ~ 60 nm. For the 



so-called "restricted primitive model" (RPM) of a electrolyte, which exhibits a 
vapor-liquid phase transition in a system of oppositely charged spherical ions of 
the same size, the Tolman's length is predicted to be twice as large. The corre- 
lation length must be sufficiently smaller than the size of the drop, allowing a 
thermodynamic approach while retaining the concept of the droplet interface. 
Therefore, in these systems, for a micron-size drop, which at AT = 10 -4 is still 

much larger than the correlation length/thickness of the interface, the curvature 
correction to the surface tension is about 1-2%. Faraway from the critical point, 
Tolman's length is only a small fraction of the molecular size and, in order to 
obtain a 1-2% correction to the surface tension, one needs a nanosize drop. 

In conclusion, let me emphasize the principal result given by Eq. (13) and 
presented in Fig. 1. Tolman's length diverges at the critical point much more 
strongly than is commonly believed. In fluids exhibiting strong asymmetry of 
phase coexistence, Tolman's length may become mesoscopic near criticality, thus 
noticeably affecting the surface tension of curved interfaces. 

Valuable discussions with Edgar M. Blokhuis, Michael E. Fisher, Michael 
R. Moldover, Jan V. Sengers, Jan Thoen, and Benjamin Widom are highly 
appreciated. I am also thankful to Heather St. Pierre for assistance. 
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Fig. 1. Tolman's length (multiplied by factor of 2), calculated from Eq. (13) 
with cs = 2/3, as a function of distance to the critical temperature. The solid 
curves represent n-heptane (£q ~ 0.11 nm, 03 ~ 0.37, b 2 — 0.09, B ~ 1.84, 
An ~ 22.6, and B cr ~ 18.4 HUES]), n-hexadecane+nitrobenzene (C16-NB) 
0.18 nm, a 3 ~ 0.48, b 2 ^ 0, and B a ~ 2.4, [22 EH), and the restricted 



(Co 



primitive model (RPM) (£ 
62 ^ -0.48, B ~ 3.64, Aq 



0.18 nm for the 0.5 nm ion diameter, a$ ~ 0.14, 
11.6, and B CI ~ 4.9 HUESj). The dotted curve 



represents the positive contribution to —25 for n-heptane, diverging as 



while the dashed curve is the negative contribution, diverging as 
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